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Theoretical considerations on multiparton interactions in QCD 
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Abstract: We investigate several ingredients for a theory of multiple hard scattering in hadron- 
hadron collisions. Issues discussed include the space-time structure of multiple interactions, their 
power behavior, spin and color correlations, interference terms, scale evolution and Sudakov loga- 
rithms. We discuss possibilities to constrain multiparton distributions by lattice calculations and by 
connecting them with generalized parton distributions. We show that the behavior of two-parton 
distributions at small interparton distances leads to problems with ultraviolet divergences and with 
double counting, which requires modification of the presently available theoretical framework. 



1 Introduction 

In hadron-hadron collisions at very high energies, sev- 
eral partons in one hadron can scatter on partons in the 
other hadron and produce particles with large trans- 
verse momentum or large invariant mass. The effects 
of such multiparton interactions average out in suffi- 
ciently inclusive observables, which can be described 
by conventional factorization formulae that involve a 
single hard scattering. However, multiple interactions 
do change the structure of the final state. They have 
been seen at the Tevatron [TJ [2] and are expected to 
be important for many analyses at LHC [3j IH [5] . 

The phenomenology of multiparton interactions re- 
lies on models that are physically intuitive but involve 
significant simplifications. A brief review of the sub- 
ject can be found in '6' and an overview of implemen- 
tations in event generators in [7^. So far a systematic 
description of multiple interactions in QCD remains 
elusive. In this letter we report on some steps to- 
wards this goal. We will see to which extent the cross 
section formulae currently used to calculate multiple- 
scattering processes can be justified in QCD and to 
which extent they need to be completed. At a more 
fundamental level, we find that there is an unsolved 
problem of double counting between single and multi- 
ple hard scattering. 

We consider the case of two hard scatters at par- 
ton level. For definiteness we analyze the production 
of two electroweak gauge bosons with large invariant 
mass (7*, Z or W) and indicate which of our results 



can be generalized to other final states such as jets. 
Since the main interest in multiparton interactions is 
driven by the need to understand details of the final 
state, we keep the transverse momenta of the produced 
gauge bosons differential, rather than integrating over 
them. For the production of a single boson there is a 
powerful theoretical description based on transverse- 
momentum dependent parton densities fSl HI [TOl [11] , 
which we aim to extend to the case of multiparton in- 
teractions. Integrating over transverse momenta gives 
the more familiar formulation in terms of collinear par- 
ton distributions. Detailed derivations of our results 
and further discussion will be given in |12) . 



2 Tree- level analysis 

We begin by sketching the derivation of the cross sec- 
tion formula for double parton scattering at tree level. 
For definiteness we take two colliding protons and con- 
sider the case where the two partons coming from one 
of the protons are quarks. The corresponding graph is 
shown in Fig. [1^, which also specifies our assignment 
of momentum variables. 

We use light-cone coordinates v"^ — {v^ ± v^)/^/2 
and V = (v^,v'^) for any four- vector v and choose a 
reference frame where p moves fast to the right and p 
fast to the left, with transverse momenta p = p = 0. 
We consider kinematics where the invariant masses of 
the bosons are large and where their transverse mo- 
menta are much smaller, i.e. we require qr Q with 
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Figure 1: Graphs for the production of two gauge bosons by double (a) or single (b) hard scattering. The 
dotted line denotes the final-state cut. The decays of the gauge bosons into fermion-antifermion pairs are not 
shown for simplicity. 



ql^ ql^ 9t and ql ^ ql ^ . 

The lower blob in Fig. [1^ is described by the corre- 
lation function for two quarks in a proton, 
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Here T denotes time ordering and T anti-time order- 
ing, as appropriate for fields in the scattering ampli- 
tude or its conjugate. The association between mo- 
menta and field positions becomes clear if one rewrites 
ziki + Z2k2-yr = {y + ^zi){ki- ^r) + ^Z2ik2 + ^r) + 
\z2{k2 — \r) — {y — ^zi){ki + ^r) in the Fourier ex- 
ponent. For now we gloss over the flavor and color 
structure of $, which will be discussed in Sections [5] 
and [6l Throughout this work we consider unpolarized 
protons, so that an average over the proton spin in (JlJ 
is understood. The correlation function $ for two an- 
tiquarks in a proton is defined in analogy to ([T]), with 
interchanged roles of the q and q fields. The contribu- 
tion of graph [T^ to the cross section then reads 



d^qi d^q2 S App 



X (27rf / d'^r d^f 5'^^\r + f) 
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where Hi(ki, ki,r, f) is the squared amplitude for each 
of the two hard-scattering processes. Repeated Dirac 
indices ai, /3i, etc. are to be summed over. The sta- 



tistical factor 5 is 2 if the final states of the two hard 
scatters are identical and 1 otherwise. 

To proceed, we make the same approximations as in 
processes with a single hard scattering: 

1. use that the minus- momenta of right- moving par- 
tons and the plus- momenta of left-moving partons 
are small compared with the large scale Q. The 
constraint 5^'^\r+f) forces r+ and f~ to be small, 
although by general scaling arguments they could 
both be large, and the constraint 5'^^\qi — h — ki) 
leads to kf ss qf and k~ ~ q~ . Defining 



=s q^ and k^ 



(3) 



we thus find that the parton momentum fractions 



Xi and (fcj ± )/p 



fixed by the final-state kinematics. Note that this 
does not hold for the transverse parton momenta, 
since q^ receives contributions from both ki ± 
and ki ± if. 

2. neglect small plus- and minus-components, as well 
as all transverse momenta in the squared parton- 
level amplitudes Hi, which then depend only on 
qf and q~ , i.e. on quantities of order Q. After 
this approximation, the integrations over k~ and 

in ([2]) only concern the factor $, and those over 
k^ and only concern the factor 

3. perform a Fierz transformation for the index pairs 
{ai,l3i) and (ai,/3i) in The correlation func- 
tion $ is then multiplied with Dirac matrices that 
can carry Lorentz indices. One has to retain only 
those terms with the maximum number of plus 
indices, since the momenta on which $ depends 
are largest in the plus-direction. For $ the domi- 
nating terms have the maximum number of minus 
indices. 
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After these steps we find that the two-quark distribu- 
tions required to describe the graph [T^ read 

Fa^^a2{xi,X2,ki,k2,r) = {2nf2p+Jdk:[dk2dr^ 



r+=0 



— 0^+ 



2p^ 



nf (P^i i(x z~p+-z k ) 



X dy-d'y e^y^ {p\Oa2iO,Z2)OaAy,^i)\p) (4) 



(5) 



with bilinear operators 

Oa{y, z) = q{y - iz) Fa q{y + \z) 



z+=y+=0 

Here a = q, Aq, Sq labels the quark polarization and 

r.^h^, rA.^h'-Ts, n,^'2'^'+To (6) 

with J = 1,2. The operators in ([SJ are well-known 
from the definitions of single-parton densities for unpo- 
larized, longitudinally polarized and transversely po- 
larized quarks, see e.g. [131 [14]. Analogous definitions 
hold for antiquarks and for a left-moving hadron. The 
cross section in © can finally be written as 



1 



dxdxd^a S ^ E^Miqf) 

O'l :0,2—q.A.q,5q 
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where here and in the following we write F{xi,ki,r) 
instead of 2:2, fei, ^2, for brevity. The parton- 

level cross sections are given by 



2qj 



with quark spin projectors Pq{k) = ifc"'"7~, P\q{k) — 
\l5k^l~i Pgqik) = and corresponding 

antiquark spin projectors Pg- It is understood that 
for each = 5q both Paj.aa and ffi^aa carry extra in- 
dices j associated with the direction of the transverse 
quark polarization. Corresponding remarks hold for 
Ri = Sq. 

The difference r of transverse parton momenta can 
be replaced by the Fourier conjugate position y, both 



in the distributions 

Pai ,a2 : ki , 2/) 



(27r) 



-^ai ,a2 ('^'i 1 i 



x2p+J dy- {p\Oa2iO,Z2)OaAy,zi)\p) (9) 
and in the cross section 



da\Ti, 
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n- 1 dxi dxi d^Qi S A X L 1 

X J d^k.d^h S^^\q, - k, - 



d^y Fa^,a2{xi,ki,y) Fa^,a.^{xi,ki,y) . (10) 



Integration of the cross section over q^ and leads to 
coUinear (i.e. transverse-momentum integrated) two- 
parton densities 

Fa,.a2{xi,y) ^ y d^fcid^fc2Pai,a2(a;<,fci,y) . (11) 

The corresponding cross section formula is the basis 
for the phenomenology of multiple interactions and has 
been used for a long time. It was derived in [15] for 
scalar partons and in [16 for quarks, in ways similar 
to the one we just sketched. 

As we will discuss in Section [13] the integral in (jlip 
diverges logarithmically at order and requires ap- 
propriate regularization. Up to this caveat, which 
also applies to single-parton densities, F{xi,y) gives 
the probability for finding two partons with momen- 
tum fractions xi and X2 at relative transverse dis- 
tance y in a proton. By contrast, F{xi, ki, y) depends 
on both transverse-momentum and transverse-position 
arguments for the quarks and does not admit a prob- 
ability interpretation due to the uncertainty principle. 
Instead, F{xi, ki, y) has the structure of a Wigner dis- 
tribution [17] for the transverse degrees of freedom: 
the integral ([TT]| gives the probability to find two par- 
tons at transverse distance y, whereas the integral 
/ d'^y F{xi,ki,y) = F{xi,ki,r = 0) gives the prob- 
ability to find two partons with transverse momenta 
fei and k2. A similar discussion for generalized par- 
ton distributions can be found in ,18^. Using ([5]) and 
(|n]) we can identify ki and k2 as the "average" trans- 
verse momenta of the two partons and y as their "av- 
erage" transverse distance, where the "average" refers 
to the scattering amplitude and its complex conjugate 
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(see Fig. [l}i). In this sense the cross section formula 
ITUl) describes two hard-scattering processes where a 
quark and an antiquark with average transverse mo- 
menta ki and ki annihilate into a gauge boson with 
transverse momentum q,j = ki + ki. The two annihila- 
tion processes occur at an average transverse distance 
y, which equals the average transverse distance of the 
two quarks in one proton and of the two antiquarks 
in the other. If the cross section is integrated over qi 
and the specification of "average" for the distance 
y can be dropped. 

It is gratifying to find such an intuitive physical 
interpretation, which we will further develop in Sec- 
tion m On the other hand, we emphasize that ([TOl) 
and its qj-integrated version can be derived from Feyn- 
man graphs using standard hard-scattering approxi- 
mations, without any need to appeal to classical or 
semi-classical arguments. 

One can also Fourier transform the two-parton dis- 
tributions w.r.t. the transverse momenta fe,, 



Fai,a2 {Xij -^i, U) 



f[ I d'k.e^^''" 



= 2p+ j dy-^^ ^^{x,z-+x,z-)p+ 
X {p\ OaA0,Z2)OaAy,Zi)\p) . 



(12) 



This form is most suitable for the resummation of Su- 
dakov logarithms, as is well-known for single-parton 
distributions '8'. Indeed, the cross section formula 
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Yli^i dxi dxi d'^Qi S 



r 2 



E 

ai,a2—q,Aq,Sq 
0.1 ,a2—q,A^q,Sq 



X Jd'^yFa,^aAxi,Zt,y)Fa^^aA^i,Zt:y) (13) 

is reminiscent of the one for single Drell-Yan produc- 
tion [19]. Up to terms of order as related to the regu- 
larization of logarithmic divergences, the collinear dis- 
tributions F{xi,y) are obtained from F{xi, Zi,y) by 
setting Zi = 0. We note that by taking the complex 
conjugate of ^ one easily finds that F{xi,ki,y) is 
real valued, whereas F{xi, ki,r) and F(xi, Zi,y) can 
have complex phases. 

The graph in Fig. [1^ involves a two-quark distri- 
bution -F'ai.a2 in proton p. The cross section receives 
further contributions where the two partons in proton 
p are both antiquarks, or where one is a quark and the 



other an antiquark. The definitions of quark- antiquark 
distributions -^01,02 f-nd ^01,02 ^-nd associated cross 
section formulae are close analogs of the expressions 
given above. There are however further contributions, 
which will be discussed in Section [5l 

The preceding results can be generalized to hard- 
scattering processes initiated by gluons, such as gluon- 
gluon fusion into a Higgs boson via a top-quark loop. 
Two-gluon distributions are defined as in ([9|) , but with 
an extra factor of l/{xip^ X2P^) on the r.h.s. and with 
operators that are bilinear in the gluon field strength 

0,{y, z) = W/G+'\y - \z) G+^y + \z) , (14) 
where a — g, Ag, dg are polarization labels and 



it' xj'i 



33 



5" , 



(15) 



The operators Og and ©Ag appear in the usual den- 
sities for unpolarized and longitudinally polarized glu- 
ons. 0\^g describes linear gluon polarization (or equiv- 
alently gluon helicity flip by two units) and has been 
discussed in [201 121] ■ 



3 Power behavior 

A pair of electroweak gauge bosons can be produced by 
two hard scatters, but also by a single one. An example 
graph is shown in figure [T]d, and the corresponding 
cross section formula reads 



da 



d^kd^k 



dxi dx^ d^q^ dxi dxi d'^qi 
X 6^^\qj^ + q2-k-k)fg{x,k)fq{x,k 



(16) 



where x = xi + X2, x = xi + X2 and a is the cross 
section for qq annihilation into two gauge bosons. 
Dimensional analysis of (O and (ITCl) reveals that 



da 



1 



(17) 



for both the single and double hard-scattering mech- 
anisms. Here the small scale A^ represents q^ or 
the scale of non-perturbative interactions, whichever 
is larger. To obtain (flT)) one uses that parton distribu- 
tions do not depend on the hard scale (except for 
logarithms, which are discarded when counting pow- 
ers), whereas hard-scattering cross sections are inde- 
pendent of A^. We thus obtain an important result: 
multiple hard scattering is not power suppressed in 
cross sections that are sufficiently differential in trans- 
verse momenta. 
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The situation changes when one integrates over Qi 
and In. the double-scattering mechanism both 

transverse boson momenta are restricted to order A 
since they result from the transverse momenta of the 
annihilating partons. For a single hard scattering one 
has + 52! A, whereas the individual transverse 
boson momenta can be as large as Q. Due to this phase 
space effect one obtains 



nt 



A! 



1 



n- 



(18) 



In the transverse-momentum integrated cross section 
multiple hard scattering is thus power suppressed. 
This is in fact required for the validity of the usual 
coUinear factorization formulae, which describe only 
the single-hard-scattering contribution. 

The power behavior just discussed remains the same 
for parton-level processes initiated by gluons instead of 
quarks, and for final states other than gauge bosons. 

4 Impact parameter 

The distributions F{xi, zi, y) depend on spatial trans- 
verse coordinates for the quarks but still refer to 
a proton with definite (zero) transverse momentum. 
A representation purely in impact parameter space 
can be obtained using the methods of [221 HSl 
where impact parameter densities for a single parton 
are constructed from generalized parton distributions. 
To this end we first define non-forward distributions 
F{xi, Zi,y; A.) as in ([T^ but between proton states 
and — ^A) with different transverse mo- 
menta. Introducing the wave packet 



-ibp \„+ 



\P ,P) , 



(19) 



which describes a proton with definite transverse po- 
sition 6, one can show that 



2p 



^ 27r 27r 




Xi 



y 



7 



X2 



Figure 2: Visualization of the cross section formula 
((22|) when and <72 are integrated over. 



see [24]. Integrating F{xi, Zi,y;b) over b one recovers 
the distributions F{xi, Zi,y), so that the cross section 
(IT3I) can be cast into the form 



1 



dx dx d^a s ^ U'^^^^.^M) 

lli^iaXiaXta q^ ^ ai,a2=q_Aq_M_'- '^='^ 
Q-i ,o,2—q,A.q.Sq 



d^Zj 



d^y d^b d^b 



-x Fa,,a2 {xt ,Zt,y;b) Fa,,a2 {xt ,Zt,y;b), 



(22) 



which has a simple geometric interpretation in impact 
parameter space. Taking the average of transverse po- 
sitions in the amplitude and its conjugate as in Sec- 
tion [21 one identifies y as the average distance between 
the two scattering partons, b as the average distance 
between parton 2 and the right-moving proton, and 
b as the average distance between parton 2 and the 
left-moving proton. This is illustrated in Fig. [2l for the 
case where the cross section is integrated over q^, so 
that Zi = and the positions in the amplitude and its 
conjugate coincide. 

The <7i -integrated version of ([2^ has previously 
been derived in [25^ . Given the work in [TSl [H] and 
[25] , we disagree with the statement in [26] that the im- 
pact parameter picture of multiparton interactions has 
until now been based on "semi-intuitive reasoning" . 



{p+,-b- ^d\Oa:,{Q,z2)Oa,{y,zi)\p+,-b+ ^d) 5 Interference contributions 



XlZl - X2Z2) Fai,a 



(20) 



with 



Fai ,a2 



t,Zi,y;b 
d^A 



ibA 



i2n) 



(21) 



The difference d in the transverse positions of the two 
proton states is a consequence of Lorentz invariance. 



The multiparton distributions discussed so far have 
an interpretation as probabilities or as pseudo- 
probabilities in the sense of Wigner distributions. 
However, there are contributions to the cross section 
which involve distributions that are interference terms 
rather than probabilities. Figure [3^ shows an exam- 
ple where the parton with momentum fraction zi is a 
quark in the scattering amplitude and an antiquark in 
the conjugate scattering amplitude. Such interference 
terms in the fermion number of the partons have no 
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So far we have not paid attention to the quark flavor 
structure of the double-scattering process. One read- 
ily finds that there are also interference terms in the 
flavor quantum numbers of the partons. An example 
is given in Fig. [5]d, where the parton with momen- 
tum fraction xi is a u quark in the amplitude and 
a d quark in the conjugate amplitude. The corre- 
sponding distribution is given by the matrix element 



Figure 3: Example graphs for interference terms in 
fermion number (a) and in quark flavor (b). The blobs 
indicating two-parton distributions are not shown for 
simplicity. Labels q and q indicate whether a line is 
represented by a quark field or a conjugate quark field 
in the corresponding matrix element. Momenta are to 
be assigned as in Fig. [T^. 

equivalent in single hard-scattering processes, where 
they are forbidden by fermion number conservation. 
Their contribution to the cross section has the same 
structure and in particular the same power behavior 
as the contributions discussed in Section [51 To the 
best of our knowledge such interference terms are not 
included in existing phenomenology. 
Let us introduce a shorthand notation 



r 2 



n 



d^z 



(27r)3 

X 2p+ J dy^ {p\ ip{y - \zi)ip{-\z2) 



(23) 



for the Fourier transformed matrix element of a prod- 
uct of field operators (p, with indices assigned as 

1 O J/ + -s-7- mom. fract. Xi in amplitude 

2 o 

3 O 



-r^ mom. fract. X2 in amplitude 

mom. fract. X2 in conjugate ampl. 
4 o y — -s-)- mom. fract. xi in conjugate ampl. 



The numbering corresponds to the order of the parton 
lines in figures [T^ and [U from left to right. We then 
respectively have 

Faua2{Xt,k„y) = {{{q3Ta2q2) {qA^aiqi))) , 

Fa,-a2{x^,k„y) = {{{q2T-a2q?.){qATa,qi))) (24) 

for a two-quark and a quark-antiquark distribution, 
whereas the quark-antiquark interference distribution 
associated with the lower part of figure [3^ is given by 



6 Color structure 

In contrast to single-parton densities, where two par- 
ton fields are always coupled to a color singlet, multi- 
parton distributions have a nontrivial color structure, 
which we now discuss. A color decomposition of two- 
quark distributions is given by 



F;jf,kk' = i{q3.M'^a2q2,k){qi,j'Yatqi,i))) 



2A^ 



1 



oTTi j-a j.a 



(26) 



where j, j', fc, k' are color indices and N is the number 
of colors. For brevity we omit the labels ai,a2 on 
F in this section. describes the case where the 
quark lines with equal momentum fractions xi or X2 
are coupled to a color singlet, whereas in ^F they form 
color triplets. Obviously, a probability interpretation 
is only possible for ^F, whereas ^F may be regarded as 
an interference term in color space. 

The color structure of Fjji^kk' can alternatively be 
parameterized by ^F and the matrix element 



Sj'k Sk'j Fjj> ^kk' — {{{q3J^a2'l2,k) iq4.k^aiqi,j))) 

1 



N 



F 



N 



'F, 



(27) 



where quark lines with different longitudinal momenta 
are coupled to color singlets. We note that (p7)) be- 
comes equal to *F in the large- A^ limit, except if one 
has \^F\ <^ \'^F\. By a Fierz transform one can express 
(|27|) in terms of matrix elements 

'F = (((g4,fcr,,g2,fc) (g3,,T„,gij))) , (28) 

where the bilinear quark operators have no uncon- 
tracted color or spinor indices. 

To illustrate that the color octet combination ^F 
need not be small let us consider a three-quark sys- 
tem, as it is done in constituent quark models. Since 
the color part of the three-quark wave function is ejki , 
the color structure of a two-quark distribution then 
reads 



Ia^,a2{x^,k^,y) = (( ((?2 Tg^ ^4) (^3 Tai <7l) )) • (25) 



Fjj',kk' fx ejki ej'k'i — Sjj' Skk' — ^jk' ^kj' 



(29) 
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where I is the color index of the spectator quark and 
therefore summed over. From this one readily obtains 
= —y/2{^F), whereas the combination in (l?7l) is 
found to be equal to — {^F). 

With a suitable assignment of color indices, decom- 
positions analogous to (j26p can be written down for 
distributions F involving one or two antiquaries instead 
of quarks, and for the interference distributions / in 
([25l) . Given the normalization of ®F chosen in ((26)) . 
color singlet and color octet distributions enter with 
equal weight 



^F^F + ^F^F)/N' 



(30) 



in the cross section if each hard scatter produces an 
electroweak gauge boson or any other color-singlet sys- 
tem. 

Two-gluon distributions have a more involved color 
structure. We first couple each of the gluon pairs {1,4} 
and {2,3} to an irreducible representation and then 
couple the two pairs to an overall color singlet. This 
gives 



-laa ,bb 



1 



(iV2 - 1)2 



^F 6°'°' S'''' 



N 



1 



bV 



jcaa jcbb' 

(31) 



where we use a shorthand notation HaG'' = 
Wj G'^^'^ G'^^'^ for the contraction of Lorentz indices. 
Each of the pairs {1, 4} and {2, 3} is coupled to a sin- 
glet, an antisymmetric and a symmetric octet in ^F, 
^F and ^F, respectively. The ellipsis in ([5T|) stands 
for terms where the pairs are in higher representations, 
which are 10, 10 and 27 for SU{i). The correspond- 
ing tensors in color space can be found in [27] , cf . also 
App. A of [28]. In hard-scattering processes that pro- 
duce color singlet states, the distributions enter as 



'^F^F^ 



ApAp 



...)/(iV2 - 1)2. (32) 



Of course there are also mixed quark-gluon distribu- 
tions. The quark lines can only couple to a color sin- 
glet or octet, which has to be matched by the gluon 
lines in order to obtain an overall singlet. A complete 
decomposition is thus given by 



jpaa 

^33' - 



1 



{{{q3,r^a,q2,,){Gl n,,G?))) 



1 



N{N^ - 1) 



^FS' 



33 



2N^ 

iV2 - 4 



(33) 




Figure 4: Example graph with a coUinear gluon Ic and 
a soft gluon Ig. 



We note that the color structure of two-parton distri- 
butions has already been discussed in PT' , using a basis 
where the parton pairs {1, 2} and {3, 4} are coupled to 
states of definite color. 



7 Wilson lines 

Our discussion so far has been concerned with tree 
graphs as in Fig. [T] At this level, our results can read- 
ily be generalized to other hard-scattering processes, 
in particular to jet production with the well-known 
subprocesses qq — >■ qq, qg — ^ qg, etc. 

A proper factorization formula in QCD must of 
course include corrections to the tree-level cross sec- 
tion, and in particular take care of additional gluon 
exchange. Detailed studies of factorization in terms of 
transverse momentum dependent parton distributions 
have been performed for single Drell-Yan production 
and its crossed counterparts, e^e~ annihilation into 
back-to-back hadrons and semi-inclusive deep inelastic 
scattering [1 [H [13 HI] . In [12] we argue that the fac- 
torization proof for Drell-Yan production can to a large 
part be extended to double hard-scattering processes 
producing uncolored states such as electroweak gauge 
bosons. We restrict ourselves to such processes from 
now on. For hard-scattering processes like jet pro- 
duction, serious problems for establishing transverse- 
momentum dependent factorization have been encoun- 
tered even for single hard scattering [53] ■ A treatment 
of the multiple-scattering case will probably have to 
wait until this situation is clarified. 

Starting from tree-level diagrams, there are two 
types of additional gluon exchange that are not power 
suppressed in the large scale and hence need to be 
taken into account systematically. The fist type con- 
cerns gluons which emerge from the subgraph repre- 
senting the partons in the right-moving proton and 
which attach to a hard-scattering subgraph, see Fig.UJ 
These gluons move fast to the right, and to avoid 
power suppression their polarization must be in the 
plus-direction. To leading-power accuracy, the effect 
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of these gluons can be represented by Wilson lines that 
appear in the operators defining parton distributions 
and make them gauge invariant. For gauge boson pair 
production, each quark or antiquark field in (jl2[) is to 
be multiplied by a Wilson line according to 



qj{z) [W{z,v)]jkqk{z) , 
Qjiz) qk{z) [W\z,v)]kj 



with 



W{z, v) — Pexp 



dX vA'^iz-Xv) e 



(34) 



, (35) 



where j and k are color indices, P denotes path order- 
ing, and the sign convention for the coupling g is such 
that the covariant derivative reads = d^^+igA^^^°-t°- . 
An analogous discussion holds for left-moving gluons 
in the left-moving proton. 

One naively expects the vector v to point in the 
minus direction for right-moving and in the plus di- 
rection for left-moving partons. This leads however to 
rapidity divergences in the parton distributions, due 
to contributions from left-moving gluons in a right- 
moving hadron and vice versa. To exclude this un- 
wanted kinematic region and to remove the divergence, 
one can take v with nonzero plus- and minus compo- 
nents [5J [TU] • This results in an additional parameter 



{2pvf 



(36) 



in the parton distributions for proton p. Their C de- 
pendence is connected with Sudakov logarithms and 
will be discussed in Section [TU An adequate choice of 
C in cross section formulae is the hard scale Q. 

The second type of unsuppressed gluon exchange is 
between the right- and left-moving partons as shown 
in Fig. m provided the gluons are soft and thus do not 
take partons far off shell. In processes with small ob- 
served transverse momenta in the final state, the effect 
of these gluons does not cancel, but it can be described 
by a so-called soft factor, which is defined in terms of 
vacuum expectation values of Wilson lines. Proper 
care needs to be taken to prevent double counting, be- 
cause the Wilson lines W{z,v) in the parton distribu- 
tions include soft gluon momenta as well [5J [TUl [TT] . 

We now turn to collinear multiparton distributions. 
It is well known that for single-parton distributions 
the rapidity divergences just mentioned cancel between 
real and virtual graphs [3^ , so that in this case the di- 
rection V can be taken exactly lightlike. The same ar- 
gument holds for two-parton distributions in the color 
singlet channel. As can be seen from pT|) and (fT2|) . the 
transverse separation Zi of quark and antiquark fields 
is zero in the collinear distribution F{xi,y). For color 



singlet distributions ^F, the Wilson lines whose color 
indices are contracted therefore combine to a Wilson 
line of finite length, going from y — ^zi to y + ^zi or 
from —\z2 to along the light cone. 

The same is, however, not true for color octet distri- 
butions . Rapidity divergences from real and virtual 
graphs do not cancel in that case, because compared 
to ^F some graphs change their color factors whereas 
others do not. One must hence keep v away from the 
light cone even in collinear octet distributions. Fur- 
thermore, the color indices of Wilson lines with equal 
transverse positions do not match, so that these Wil- 
son lines cannot be combined to a line of finite length. 
As a consequence, terms with color octet distributions 
in collinear factorization formulae will have a very dif- 
ferent structure from the usual one. 

It should be possible to extend the previous discus- 
sion to gluon distributions. For the Wilson lines in 
single-gluon densities a detailed analysis can be found 
in [21] , whereas to our knowledge the soft-gluon sector 
has not been elaborated. 



8 Parton spin correlations 

The cross section formulae in Section [2] have a non- 
trivial dependence on parton spin, even though they 
are for unpolarized colliding protons. Fq^q and its 
counterparts for antiquarks describe unpolarized par- 
tons, whereas all other distributions F^q^Aq, PqAq 
etc. describe correlations between the polarization of 
the two partons, or between their polarization and 
the transverse vectors y and ki. Such spin correla- 
tions have been pointed out in 27 , but to our knowl- 
edge they have not been implemented in phenomenol- 
ogy. We find that 16 scalar or pseudoscalar func- 
tions are required to parameterize the spin structure of 
Fai.a2i^ijf^iiy) a givcu quark flavor combination. 
Eight scalar functions remain when one integrates over 
ki and k2 to obtain i^Qj^ Q^ (xi, y). 

To see that parton spin correlations in the proton 
need not be small, let us consider a SU{6) symmetric 
three-quark wave function. As is well known, this gives 
Au/u = 2/3 and Ad/d = — 1/3 for the average longitu- 
dinal polarization of u and d quarks, which reproduces 
the trend observed in polarized quark densities at x 
around 0.3. Similarly, one obtains Fau.Au/ Fu,u = 1/3 
and FAu,Adl Fu,d = —2/3 and thus a considerable cor- 
relation between the longitudinal polarization of two 
quarks. 

Parton momentum fractions in processes at LHC are 
typically well below 0.3, say of order 10~^ or smaller. 
In that region polarized single-parton distributions (to 
the extent that they are known) are small compared 
with their unpolarized counterparts. This means that 
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there is httle correlation between the respective polar- 
izations of a parton and the proton when they are far 
apart in phase space. From this we should however 
not conclude that distributions like -FAg,Ag are unim- 
portant at small but comparable Xi ^ X2, since they 
describe a correlation between two partons that are 
relatively close in phase space. 

Parton spin correlations have important conse- 
quences in multiple interaction processes. Let us again 
consider gauge boson pair production. For longitu- 
dinal polarization of both the quark and the anti- 
quark that annihilate into a boson, each spin pro- 
jector FAg and Fa^ in the hard-scattering cross sec- 
tion ([5]) contains a 75. After anticommuting one 75 
through the fermion trace, one obtains a factor 7I = 1, 
so that the cross section depends on the combination 
Fq,q Fq,q + F^q^Aq ^A(j,Ag of two-partou distributlous. 
Longitudinal spin correlations between quarks and an- 
tiquarks thus directly affect the rate of gauge boson 

pair production by multiple scattering. 

■ -/ 

The distribution Fg^ g^ gives rise to even more strik- 
ing effects. Its decomposition into scalar functions con- 
tains a term with S^^ , which describes a correlation 
proportional to the product S1S2 of the two transverse 
quark polarization vectors. Recall that the quark field 
bilinears qia'^^j^q defining the distribution Fg^ g^ are 
chiral odd, which corresponds to quarks with opposite 
helicities in the amplitude and its conjugate. As a re- 
sult, transverse quark polarization does not contribute 
to the production of W bosons. There is however a 
term with Fg^ g^ Fg^g^ in the cross section for produc- 
ing two neutral gauge bosons, 7*7*, 7*^ or ZZ. This 
term gives in particular rise to an angular correlation 
proportional to cos 2ip, where ip is the azimuthal angle 
between the momenta of the negatively charged lep- 
tons or the quarks from the boson decays. We thus 
obtain the important result that correlations between 
transverse quark and antiquark polarizations give rise 
to a correlation between the decay planes of the two 
produced bosons. Moreover, this type of correlation 
does not arise if the two bosons are produced by a 
single hard scattering such as in Fig. [T|d, where for 
Qt Q one obtains an angular modulation with cos 
but none with cos 2Lp. 

It is natural to expect that spin correlations between 
partons in the proton also result in correlations be- 
tween the scattering planes in other double-scattering 
processes, such as the production of two dijets. We 
note that in [32 the absence of such correlations was 
explicitly assumed, which in view of our discussion may 
not be adequate. It would also be interesting to study 
whether polarization induced angular correlations can 
contribute to the so-called "ridge effect" observed in 
pp collisions by CMS [33] . 



9 Mellin moments 

If one takes Mellin moments / dx^ a;"^ / dx2 (^F) 
of color singlet distributions, the light-cone operators 
Oa{y, z) in ([5]) turn into local twist- two operators that 
are well known from the operator product expansion. 
(As follows from the discussion in Section [71 the cor- 
responding operators for color octet distributions still 
contain Wilson lines and will not be discussed here.) 

Let us consider two unpolarized quarks and take the 
lowest Mellin moment in both xi and X2, 



Mq,q{y^) 



1 /•! r 

dXi / dX2 ^Fq,q{xi,X2,y) 

Jo L 

^Fg^g{Xl,X2,y) - '^Fg^g{Xl,X2,y) 



+ ^Fg^q{xi,X2,y) 



p^ 



dy-{p\Og{0,0)Og{y,0)\p). (37) 



We introduce the local operator 0^{y) = q{y)^^q{y) 
and decompose 

{p\0''{Q)0^{y)\p) ^2p^p^ {00){py,y^) + ■ ■ ■ , (38) 

where the ellipsis represents terms with uncontracted 
vectors j/^, y'^ or the metric tensor g^'^ . The reduced 
matrix element (CO) can only depend on the invari- 
ants py and y^ . We now choose a frame where p = 



and y^ = 0, so that py 



p+y- 



and y — —y . Using 



Oq{y,0) = ^0^{y) \y+=Q we can then write ([37)) in a 
manifestly covariant form 



Mq.,q{y^)= / d{py) {00){py,y^) 



(39) 



It is straightforward to generalize this procedure to 
higher Mellin moments and to the other quark polar- 
izations Aq and 6q. In each case one can express the 
Mellin moment in terms of the matrix element of a 
product of two local twist-two operators. 

The matrix element in p9|) can be evaluated on a 
lattice in Euclidean spacetime if one takes = 0. This 
is rather similar to recent lattice studies of transverse- 
momentum dependent single-quark distributions [341 
[35]. The restriction to y^ — entails 



(pyf/i-y^) = {pyf/f < p' 



(40) 



where p and y denote the spacelike three- vectors. Re- 
sults from a discrete Euclidean lattice are hence not 
sufficient for evaluating the integral over all py at fixed 
y^ in p9|) . This is completely analogous to the situa- 
tion discussed in [35j . Despite this limitation of prin- 
ciple, we hope that lattice data in a certain range of 
py and y^ will one day provide genuinely nonpertur- 
bative information about the behavior of multiparton 
distributions. 
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10 Connection with generalized 
parton distributions 

To develop a viable phenomenology of multiple inter- 
actions, one needs a simple ansatz for multiparton dis- 
tributions that can be progressively refined. For the 
distributions which admit a probability interpretation, 
a natural starting point is to replace them by the prod- 
uct of single-parton densities. 

To formalize and extend this ansatz, we insert a 
complete set of intermediate states |X) (X| between the 
operators Oaa and Oa^ in the definition of two- 
parton distributions in impact parameter space. This 
gives a product of single-parton operators sandwiched 
between a proton state and X . If we assume that the 
ground state dominates in the sum over all X and take 
the intermediate proton states in the impact parame- 
ter representation then F{xi, Zi,y;b) involves a 
product of single-proton matrix elements 

(p+,-6~id|a,(0,Z2)|p'+,b') 
X (p'+,6'|a,(2/,zi)b+,-6+id) 

X (p'+,6'-y|a,(0,zi)|p+,-b-y-t-id). (41) 



which has a straightforward physical interpretation as 
sketched in Fig. O In different guises, this relation is 
at the basis of most phenomenological studies and has 
long been used in the literature, see e.g. [36l [37l [38] 
and [221 HD] ■ As was noticed in Pl] , the convolution in 
(P5|) simplifies to a product if one Fourier transforms 
to transverse momentum space, where one has 



(46) 



with fa{x; A) = / d^b e*^/a(x; 6). 

In the previous argument we have neglected the spin 
of the proton. When inserting intermediate proton 
states between and Oa^ in a two-parton distribu- 
tion for an unpolarized proton, we schematically have 

i ^ (p, A I a, I p', A') (p', A' I a, I P, A) (47) 



A 

X 



A, A' 



with proton helicities A and A'. The sum on the r.h.s. 
includes matrix elements where the proton helicity dif- 
fers in the bra and the ket state. As an example let us 
consider the coUinear distribution for two unpolarized 
quarks. We then find 



Integrating the phase factor (p'-p)^ over y- sets ^Fq,q(x,,r) w ff«(x2, 0, -r^) iJ'^(xi, 0, -r^) 
p'+ — p+ in ((20)) . and we obtain a representation 



^Fai,a2{xi,z,,y;b) « /q^ (x2 , 2:2 ; 6 -f- ^XiZi) 

^ fat{xi,zi;b + y - \x2Z2) (42) 

in terms of single-quark distributions 



Ami 



£;«(x2,0,-r2)^9(xi,0,-r2), (48) 



fa{x,z;b) 



dz 



27r 

X (p+,iA|a(0,z)|p+,-iA) (43) 
in impact parameter space. These distributions satisfy 

^ e"^-p+ _6 - \d |a(0, z) -b+\d) 

(44) 



^6^^\d- xz) fa{x,z;b) 



in analogy to ((20)) . For z — they are the impact pa- 
rameter dependent parton densities fa{x;b) discussed 
in [22J [23] and give the probability to find a parton 
with momentum fraction x and distance b from the 
proton center. 

At Zi = the relation ((42)) involves only collinear 
distributions. Integrating over b we get 

'^Fa,,a2{xi,y) ~ J d^b fa2{x2;b) fa,{xi;b + y) , (45) 



where the first term corresponds to A = A' and the sec- 
ond one to A = —A' in ((TT)) . The generalized parton 
distributions (GPDs) H''{x,^,t) and Ei{x,^,t) can be 
studied in hard exclusive processes such as deeply vir- 
tual Compton scattering or vector meson production. 
For their definitions and further information we refer 

to gusg. 

We emphasize that the relations ([42]), (gSJ and ((46)) 
are obtained by restricting a sum over all intermediate 
states to a single proton in a selected helicity state. We 
do not have a justification or a strong physical motiva- 
tion for this restriction, other than stating a posteriori 
that it is tantamount to neglecting any correlation be- 
tween two partons in the proton. It seems plausible 
to assume that this is a reasonable first approxima- 
tion, at least in a certain region of variables, but one 
should not expect such an approximation to be very 
precise. The result ((48)) already goes beyond the com- 
plete neglect of correlations, and the relative size of 
the term depending on may be taken as an indica- 
tor for the level of precision of the simplest factorized 
ansatz. Possible deviations from ((45)) and their conse- 
quences for multiple scattering cross sections have e.g. 
been discussed in [40ll43ll44] . 
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Figure 5: Illustration of the approximation (j45l) of a two-parton distribution in terms of single-parton distribu- 
tions. Implicit in the figure is the representation of these distributions as squares of light-cone wave functions 
in the proton. 



So far we have discussed the distributions 01,02: 
which can be interpreted as probabilities or pseudo- 
probabilities. The formal derivation we have sketched 
can however be extended to distributions that have 
interference character. Repeating the above steps for 
the distributions appearing in Fig.[3}D, one obtains ma- 
trix elements of operators dVu or uTd that transfer 
isospin. For a proton target, the ground state among 
the intermediate states inserted between the two oper- 
ators is then a neutron. Isospin symmetry relates the 
resulting matrix elements as {n\dTu\p) = {p\uTd\n) = 
{p\uTu\p) — {p\dTd\p). Similar relations for strange 
quarks assume flavor SU(3) symmetry and can be 
found in [411142]. 

Inserting physical intermediate states is useful be- 
tween the color singlet operators appearing in but 
not between the color octet operators qTt°^q in ^F. One 
can however repeat the preceding construction for the 
distributions ^F defined in (pS)) . In that case the two 
partons represented by a color singlet operator carry 
different plus-momentum fractions xi and X2. This 
implies p'"*" ^ p^ in the resulting product of proton 
matrix elements, which are therefore associated with 
GPDs at nonzero skewness ± (p — p')'^ / {p + p')^ ■ 
Likewise, the rearrangement of fields to color singlet 
operators in the interference distribution (|25|) leads to 
GPDs with nonzero ^. We note that in the exclusive 
processes where GPDs can be measured, one always 
has ^ 7^ because the scattered proton must have a 
smaller momentum than the proton target. 

GPDs give rather direct information about the dis- 
tribution of partons in impact parameter b, which is 
Fourier conjugate to a measurable transverse momen- 
tum A in physical processes. By contrast, the inter- 
parton distance y in multiple interactions is integrated 
over in the cross section formula (fTO|) and not directly 
related to any measurable kinematic distribution. Al- 
though the connection between GPDs and multiparton 
distributions is not exact, it provides an opportunity 
to obtain some quantitative information about impact 
parameter dependence from an independent source. 



An important example is the correlation between im- 
pact parameter and longitudinal momentum of par- 
tons, for which there are clear indications in GPDs (see 
Section 4.4.5 of [IS] and references therein) and which 
has recently been studied for multiparton interactions 
in |46]. 

11 Ladder graphs 

The factorization formulae in Section [2] are for kine- 
matics where the transverse boson momenta are much 
smaller than the large scale Q. This includes but is 
not limited to the region where is comparable to 
a hadronic scale A. In this and the following section 
we consider the region of intermediate transverse mo- 
menta A <C (7t ^ Q- If 9i and ^3 are large com- 
pared with A then at least some of the transverse mo- 
menta in the two-parton distributions must be large as 
well. This opens the way for a perturbative description 
where the transverse-momentum dependent distribu- 
tions factorize into collinear distributions and a hard 
parton-level scattering. This increases the predictive 
power of the theory, given that collinear two-parton 
distributions depend on fewer variables and are less 
numerous. Note that the scale qT of the hard scatter- 
ing in the two-parton distributions is still much smaller 
than the scale Q of the multiple hard scattering in the 
proton-proton collision. 

The type of graphs that describe the hard scatter- 
ing in F{xi,ki,r) depends on the relative size of the 
transverse-momentum arguments. For instance, the 
single- ladder graph in Fig. [6] gives large fci, whereas 
a double-ladder graph with an additional gluon ex- 
changed between the quarks with momentum fraction 
X2 leads to large fci and fc2- In both cases no hard 
parton is exchanged between the lines with momen- 
tum fraction xi and those with momentum fraction 
X2, so that r is of order A. This corresponds to an 
interparton distance y of hadronic size. 

The computation of the graph in Fig. [6] proceeds 
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A-1 



/C2 + k2- I'r fci + ir 




Figure 6: Ladder graph for a two-parton distribution 
at large ki and small fc2 and r. 

exactly as for single-quark distributions at large trans- 
verse momentum and can e.g. be found in 1471 148] . 
One obtains 



J tfk2 F-^{xi,X2,ki,k2,y) 
J.J' (^\ 

\UiJ 



E 



dui 

Ul 



Fig.m 



F-'\ui,x2,y), (49) 



where P has a logarithmic dependence on k^ which we 
have not displayed for brevity. We have integrated the 
distribution over k2 since this turns out to be the quan- 
tity that appears in the cross section at large and 
^2- The indices J and J' indicate that one has a ma- 
trix structure if transitions between gluon and quark 
distributions are permitted by the quantum numbers. 
Quark-gluon transitions occur for the combinations 



SpJ ^ 




ApJ 



(50) 



where the second parton index a can indicate an unpo- 
larized or polarized quark or antiquark. If a indicates a 
gluon, one should replace ^Fg^a + ^Fq^a by ^Fq^a + ^Fq^a 
and ^Fq^a - ^Fq^a by ^Fq^a - ^Fq^a- Taking into account 
the color factors in the definitions (pS)) . (I5T1) and 
we find splitting matrices 



ipJJ' ^ 



SpJJ' _ 



ApJJ' 



npPq, 
NPnn 



P 

2N qq 



8 



Ar2 




--^ P 

2N ^qq 



2{N' 



TTrry npPqg 



8 



Pn, 



N 



(51) 



_ p 

2 ^99 



for N colors and np quark flavors, where the func- 
tions Pqq{z), Pqg(z), etc. are defined without color fac- 
tors and differ from the usual DGLAP splitting func- 
tions at most by terms proportional to i5(l — z). The 
color factors in ([5T|) agree with those in if one re- 
stores a missing factor y/2/N in the expression of Pg/ , 
eq. (54b) of that paper. This factor is also required 
for the normalization of Pgf as a projector in eq. (53) 
of [49]. 

We see that color factors are always smaller in the 
octet channels than in the singlet channel, with the 
biggest suppression occurring for Pqq. In the large- A'^ 
limit the singlet matrix ^P has eigenvalues NPgg and 
^Pqq, whereas for both ^P and "^P one of the eigenval- 
ues is ^Pgg and the other is of order 1. This suggests a 
relative suppression of color octet distributions at large 
transverse momentum, but quantitative estimates are 
needed to assess how important this suppression is. 

Quark-gluon transitions do not occur for the combi- 
nations J2q i^Fq.a — ^Fq^q], for the difference of distri- 
butions for two quark flavors, and for matrix elements 
where the quark flavors differ on the two sides of the 
final-state cut as in Fig. [Jb. In these cases the in- 
dices J and J' in (|49)) take only one value, with the 
splitting function P being CpPqq for color singlet and 
^TN^qq '^'-"^ color octet combinations. Furthermore, 
there is no transition between gluons and the inter- 
ference distributions in Fig. [5^. From the experience 
with single-parton distributions one can expect that 
at low xi and high ki those combinations of quark 
and antiquark distributions will dominate that receive 
a contribution from gluons on the r.h.s. of pQ]). 

The splitting matrices for longitudinally polarized 
quarks and gluons contain different splitting functions 
APgg, APqg, ctc. but havc the same color structure 



as (|51|). No transitions occur between 5q^ 6q and Sg, 
since the former describe parton helicity flip by one 
unit and latter by two units. A relation similar to (|49l) 
can be derived for double-ladder graphs and involves 
the same splitting matrices as above. 

We finally note that for qr ^ A the contributions 
from ladder graphs lead to a power behavior 



s da 



Yli^i dxi dxi d^qi 



ladders 



A^ 

4 



of the double-scattering cross section scaled with s^. 



12 Parton splitting 

We now consider the case where \r\ ^ \ki\ is large, 
which corresponds to perturbatively small distances 
\y\. Power counting analysis shows that the dominant 
contributions to the double-scattering cross section are 
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Figure 7: Lowest-order graph for a quark-antiquark 
distribution at large \r\ ~ \ki\. The two partons 
with momentum fractions xi and X2 originate from 
the splitting of a single gluon. 

from graphs where two partons in the t channel split 
into four partons. A lowest-order example for i^ai,a2 
is shown in Fig. [T) At next order in graphs appear 
where the hard scattering if fully connected, e.g. with 
an additional gluon exchanged between the left- and 
right-hand sides of Fig. [T] The power behavior of the 
scaled cross section is found to be 



s da 



Hi^i dxi dxi d^Qi 



splitting 



1 



(53) 



if in at least one of the colliding protons the two-parton 
distribution has a fully connected hard scattering. If 
in both distributions the hard-scattering subprocess 
is disconnected as in Fig. [HI the scaled cross section 
behaves like 1/A^ instead, but the final-state phase 
space is then restricted to \qi + q2\^ ^■ 

We see that the contribution ((52|) from ladder graphs 
at large y is suppressed by A^/g|, compared with the 
splitting contribution ([5^ at small y. This suppres- 
sion may however be compensated by a stronger in- 
crease of ([5^ when xi and X2 become small. It is 
known that in appropriate quantum number channels 
the iteration of ladder graphs leads to a growth of 
parton densities at small momentum fractions. The 
two-parton distribution at the bottom of Fig. [6] con- 
tains one ladder between the lines with momentum 
fraction xi and another between the lines with mo- 
mentum fraction X2 , whereas at the bottom of Fig. [7] 
there is a single-parton density containing only one 
ladder. To establish whether the A^/q|, suppression or 
the small-x enhancement of the ladder-graph contribu- 
tions is more important in given kinematics requires a 
detailed study. 

The calculation of the graph in Fig. [7] reveals a num- 
ber of important issues. We give here only the contri- 
bution from the transverse-momentum dependent un- 
polarized gluon distribution fg{x,k) and omit terms 
depending on the gluon Boer-Mulders function [50], 
which describes the correlation between the transverse 
momentum and the linear polarization of a gluon. In 



Figure 8: Graph for the cross section where both two- 
parton distributions -^121,52 and -Fsi,a2 (indicated bj 
boxes) involve the splitting of one into two partons. 



the color singlet channel we find 



-^ai ,a2 {'^i : i ^) 



_ as_ fgjxi +X2,ki +k2) 
Fig.E] 47r2 xi + X2 



X r. 



Xl 



Xi + X2 



{k 



k 



(54) 



with k ~ \ {ki — ^2)- For the color octet distributions 
^Fai,a.2 we obtain the same expression times a suppres- 
sion factor —l/\/N'^ — 1. We have 



-t; 



Aq.Aq 
nil' 



niss'i^)^ -<5^^'<5"'2u(l-u), (55) 



5q^5q 



whereas Ta-^^a^ — for the other polarization combina- 
tions. We see that several nontrivial and large spin cor- 
relations are generated by perturbative splitting. To 
which extent these correlations persist at nonpertur- 
bative values of y is an interesting question from the 
point of view of hadron structure. 

With the analog of ([7]) for the product i^ai,a2-P'ai,a2 
one finds that the contribution of Fig. [5] to the cross 
section for two-boson production is proportional to 



d'k, 



d^k^ 



ki{k_ -qf 



(56) 



with q = \{qi — 92) ^.nd k± = k ± ^r. Each inte- 
gral is infrared finite but has a logarithmic divergence 
at large k±. To assess the meaning of these diver- 
gences, we first recall that the derivation of the cross 
section formula ([7]) and its analogs for other channels 
assumes transverse parton momenta I fci ±ir I <C Q- As 
it turns out, the integrand in (I56p does not decrease 
fast enough to suppress the region of large k±, where 
the approximations leading to ^ are invalid. On the 
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other hand, we note that Fig. [5] can also be read as 
a graph for producing two gauge bosons Vi and V2 in 
a single hard-scattering process. For the parton level 
amplitude it represents a box graph, which indeed di- 
verges logarithmically in the ultraviolet. The full am- 
plitude for gg — >■ Vi V2 is however finite, because the 
divergences of all contributing box graphs cancel each 
other. This was noted in 51 and has long been known 
for the related process 77 — ?> 77 with on-shell photons 
[52] . A detailed analysis of the infrared behavior of 
gg — Vi V2 is given in , and explicit expressions for 
the gg ^ Z Z amplitude can be found in |53j . 

Since the definition of two-parton distributions we 
have used so far includes the splitting contribution 
([54]) . this definition is not appropriate for the cross 
section formula ([7]) , and one or both of them needs to 
be modified. It remains a task for future work to devise 
a consistent formulation that is also suitable for practi- 
cal computations. Such a formulation must obviously 
avoid the divergent integrals in (j56p . In addition it 
must address the problem of double counting in Fig. |5] 
when the contributions from single and double hard 
scattering are added in the cross section. The analo- 
gous double counting problem for multijet production 
was pointed out in ^54j . 

13 Collinear distributions 

We have encountered collinear two-parton distribu- 
tions in two different contexts, firstly in the cross sec- 
tion for multiple scattering when final-state momenta 
are integrated over, and secondly in the calculation 
of two-parton distributions at high transverse momen- 
tum. In this section we discuss the scale evolution of 
collinear distributions at leading order in as- We fo- 
cus on the color singlet sector, which is most similar 
to the case of single-parton densities. (As noted at 
the end of Section [3 collinear color octet distributions 
involve further complications related with rapidity di- 
vergences.) 

Integrating F{xi, ki, r) over the transverse momenta 
ki gives logarithmic divergences, which is evident from 
the factor l/k\ in (|49p . The definition of collinear dis- 
tributions requires subtraction of these divergences, in 
full analogy to the case of single-parton distributions. 
We now discuss the contribution of the splitting graph 
in Fig. [7] and consider unpolarized quarks for definite- 
ness. The expression ([54| behaves like 1/fc^ for large 
k and thus gives a logarithmically divergent integral. 
This divergence needs to be subtracted as well. Notice 
that the four quark lines in Fig. [7] are far off-shell if 
one of the transverse momenta r and k is large. This 
implies that (f54| describes not only the large r behav- 
ior of -F'ai,a2 fei) f) but also its behavior for large k 



at small r. 

Calculating the graphs in Figs. [6] and [7] in 4 — 2e di- 
mensions, performing MS subtraction of the ultravio- 
let divergent terms, and adding the contributions from 
self-energy graphs, one obtains the evolution equation 

^Fg,g(xi,X2,r) 



d log /i' 



bi—q,g ^ 

^^""^ du2 „ ( X2 \ I , , 

Pq.b2\ — ] Fq^b2[Xl,U2,r} 

X2 ""2 \U2/ 



E 

f>2=9,9 
1 



Pn 



X1+X2 \Xl + X2 ' 



fg{Xl + X2) 



(57) 



which is valid for any value of r. Here Pbi,b2 denotes 
the usual DGLAP splitting functions (including color 
factors). Remarkably, the y dependent distributions 
F{xi,y) evolve differently. Working in 4 — 2e dimen- 
sions, one finds that for nonzero y 



.2-2e^^2-2efc _^ (fc + ir) (fe - ir) 



(27r) 



2-2e 



(58) 



has the finite value 1/y^ at e = and hence requires 
no ultraviolet subtraction. One thus obtains a homo- 
geneous evolution equation 

^Fq^q{xi,X2,y) 



dlog /i^ 



E 

bi=q,g 



^ ^'^ dui f Xl \ , 

PqM[ — j Fbi,q[ui,x2,y) 



+ E 

b2=q,g -"^^ 



du 



— Pq,b2 \ — ) ^^9,62 (2^1, W2, y) • 
U2 \U2' 



(59) 



This result is very natural if one considers that 
^01,02 (2;^, y) is defined in terms of the product 

[q{\z^)^+q[-\z2)\ [q{y-\zi)l+q{y+\zi)\ (60) 

of two operators with light-like field separations zf — 
z\ — 0, where we have omitted the Wilson lines be- 
tween the fields for brevity. Each of the two opera- 
tors needs to be renormalized in the same way as in 
a single-parton density, but as long as the operators 
are taken at a finite spacelike distance = — no 
further short-distance singularities appear. 
It has long been known that 



Fq,q{Xi^ — Fqq(^Xi^V 0) 



/ d'^y Fq,q{xi,y) (61) 

J rcg 



14 



evolves as in ([ST]), see [S51[5S1[57] and the recent studies 
[551 [55] . In position space representation, the inhomo- 
geneous term now appears because the behavior 



(62) 



at smaU y induces a logarithmic divergence in the inte- 
gral over y, which needs to be regularized as indicated 
by the subscript "reg" in (|6Tt . We note that F{xi) 
does not directly appear in double-scattering cross sec- 
tions but is often used as an intermediate quantity for 
modeling the distributions F{xi,y). 

We have already seen in ([55)) that the contribution of 
splitting graphs to the multiparton distributions leads 
to logarithmic divergences in the differential cross sec- 
tion. An even worse divergence appears in the cross 
section integrated over Qi and which according to 
(O and (flUl) involves an integral 



cfyF{x^,y)F{x^,y) = 



^Fix.,r)F(,. 



(63) 



This is linearly divergent in y^ according to (15^ . It 
also diverges linearly in , because for large r one has 
F{xi,r) ~ log(r^//x^). The extra ultraviolet subtrac- 
tion included in the definition of F{xi, r) is insufficient 
to regulate this divergence in the cross section. In the 
recent paper |60j a finite result is obtained by imposing 
a cutoff < min(q'f,q'|) in ([S5)) . 

In this context we wish to comment on an ansatz 
often made in phenomenological studies, where the y 
dependent two-parton distributions are written as 



F{x^,y;n) ^ f{y)F{x^;^J.) 



(64) 



with a smooth function /(y) satisfying the normaliza- 
tion condition J cPy f{y) = 1. A typical choice for 
f{y) is for instance a Gaussian. This type of ansatz is 
obviously inconsistent if F{xi, y; fi) is defined from the 
product (150]) of twist-two operators, since the /x depen- 
dence of the l.h.s. is then given by the homogeneous 
evolution equation ([59)) whereas the ^ dependence of 
the r.h.s. is governed by an inhomogeneous evolution 
equation as in ()57|) . If one instead defines the y depen- 
dent distribution as the Fourier transform of F(xi,r), 
then the ansatz ([M)) is consistent regarding evolution 
since F{xi,r) evolves according to ([57)) . However, we 
do not think that is a satisfactory definition, since it 
does not cure the divergence of the integral ([55)) in 
double-scattering cross sections. An ansatz like ([M)) 
with a smooth function f{y) does not have this prob- 
lem and may be regarded as modeling a y distribu- 
tion in which the perturbative splitting contribution 
giving rise to the singularity has been removed. 

Since the ansatz is ad hoc, one cannot say what the 



correct evolution equation to be used on both sides 
of ([64| actually is. Our discussion suggests that the 
homogeneous form ()59|) may be more appropriate, at 
least for values y of typical hadronic size, which are of 
course most important when the ansatz is used in the 
factorization formula. A theoretically sound solution 
remains a task for future work. 



14 Sudakov logarithms 

As is well known, transverse momenta which are 
much smaller than the hard scale Q of a process give 
rise to Sudakov logarithms in the cross section. These 
logarithms must be resummed to all orders in pertur- 
bation theory, which for single gauge boson production 
can be done using the CoUins-Soper-Sterman formal- 
ism [19]. We extend this formalism to gauge boson 
pair production in |12] and sketch the main results of 
our analysis in the following. 

The dependence of a two-quark distribution on the 
rapidity parameter ( defined in ([55| is governed by the 
differential equation 

d&c Of) = [g(^iC,m) + g(.t2C,m) + ^(^i,m) 



(65) 



where ^F and ^F depend on Xi, Zi, y and They also 
depend on a renormalization scale /i, but we need not 
discuss this dependence here. The kernels G and K in 
(155)) already appear in the CoUins-Soper equation |S] 
for single-quark distributions. 



df{x,z;C) 
dlogC 



[GixC,ii)+Kiz,ii)]f{x,z;0. (66) 



The matrix M mixes color singlet and color octet dis- 
tributions and is /i independent, whereas the /i depen- 
dence of G and K is given by a renormalization group 
equation 



7if(a4M)) = - 



dK{z,fi) _ dGixCfJ-) 



dlog/x 



dlogn 



(67) 



and thus cancels in G + K. Both K and M are due 
to soft gluon exchange and can be defined as vacuum 
matrix elements of Wilson line operators, similar to 
those discussed in Section [7] They can only be cal- 
culated perturbatively if the transverse distances on 
which they depend are sufficiently small. 
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The general solution of (|55l) can be written as 



with a power 



'S(xiC,Zi,Z2)-S{x2C.Zi,Z2) 



X e 
with 



LM(zi,Z2,y) 



K{zi,iMo) + K{z2,Ho) , xC 

^ log — 

2 Mo 



Mo 



M 



lK{as{fj)) log G(m,m) 

M 



and 



L = log 



Mo 



(69) 



(70) 



The scale /xq specifies the initial condition of the dif- 
ferential equation (j65p . with a natural choice being 

/Lto OC 1/ v/Nl|N2|. 

The leading double logarithms of C/mo in come 
from the second line in (j69p . whereas terms involving 
X and M only contain single logarithms. The Sudakov 
exponent S also appears in the solution of (|66p for 
single-quark distributions [8], 

/(a;,;2;C)=e-^("'^^"^"' ^"(x,;^). 

We thus obtain the important result that to double log- 
arithmic accuracy the Sudakov factor for a multipar- 
ton distribution is the product of the Sudakov factors 
for single parton densities, both for color singlet and 
color octet distributions. A non-trivial cross talk be- 
tween all partons, and in particular a mixing between 
color singlet and octet distributions occurs however at 
the level of single logarithms, which are known to be 
important for phenomenology. When the parton dis- 
tributions F(xi, Zi,y\ Q are inserted into the cross sec- 
tion formula ([TB]). the typical size of j^ij is l/|qi| and C 
should be taken of order Q, so that logarithms of C/ /^o 
turn into logarithms of Q /qx- 

If all distances Zi and y are small, one can calcu- 
late K and M in perturbation theory. We give ex- 
plicit results in |12j and only mention some of their 
features here. The off-diagonal elements in the ma- 
trix turn out to be color suppressed, but only 
by In the limit where \zi\ <^ \y\ we find that 
the off-diagonal elements are additionally suppressed 
by \zi\\z2\/y'^ and that the diagonal element for the 
color octet is smaller than the one for the color singlet. 
This results is a suppression of octet distributions by 



Sf(x„z„y;C) 
iF(a;,,2;,,y;C) 



l^l||^2| 

y2 



A = mm 1, A* — log 

TT ^0 



(73) 



(68) In the cross section one has \zi\ ~ l/l^il, so that for 
large qt the factor ([7^ disfavors color octet distribu- 
tions in a wide range of y. To study quantitatively 
the importance of color-octet suppression, one needs 
to extend the perturbative result ([72]) to the region 
of non-perturbative distances y. This has not been 
achieved yet. 

15 Conclusions 

We have studied several aspects of multiparton in- 
teractions in hadron-hadron collisions. Our theoret- 
ical framework is hard-scattering factorization, which 
requires a large virtuality or momentum transfer in 
each partonic scattering process but is valid in the full 
range of parton momentum fractions. A complemen- 
tary approach, based on the high-energy limit and us- 
ing BFKL methods is discussed in [6ll |62l [63] . 

The basic cross section formula for multiple inter- 
actions can be derived at tree level using standard 
hard-scattering approximations and has an intuitive 

(71) geometrical interpretation in impact parameter space. 
We have shown that it can be formulated at the level 
of transverse-momentum dependent multiparton dis- 
tributions, which permits a description of the trans- 
verse momenta of the particles produced in the hard 
scattering. This is particularly important because it 
is in transverse-momentum dependent cross sections 
that multiparton interactions are not power suppressed 
compared with single hard scattering. 

The derivation of the tree-level formula for double 
hard scattering exhibits nontrivial contributions from 
correlation and interference effects. They have partly 
been discussed earlier in the literature |27] but are not 
included in current phenomenological models. The po- 
larizations of two partons can be correlated even in an 
unpolarized proton, and we have shown that such cor- 
relations can significantly affect both the overall rate 
and the angular distribution of final-state particles in 
multiple interaction processes. Two-parton distribu- 
tions have a nontrivial color structure since each par- 
ton can carry a different color in the scattering am- 
plitude and its complex conjugate. In addition, there 
are interference terms in fermion number and flavor as 
shown in Fig. [3l 

To develop a reliable phenomenology, one needs in- 
formation about the size and kinematic dependences of 

(72) two-parton distributions. One can relate them to gen- 
eralized parton distributions for single partons, which 
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are experimentally accessible in exclusive scattering 
processes, but this requires an approximation whose 
reliability we cannot quantify. Nevertheless this re- 
lation offers some guidance, especially about the in- 
terplay between longitudinal momentum and trans- 
verse position of partons, as well as the size of two- 
parton distributions that are of interference nature and 
hence do not have an intuitive probability interpre- 
tation. For transverse parton momenta above a few 
GeV one can compute transverse-momentum depen- 
dent distributions in terms of collinear ones and thus 
has more predictive power from theory. We find a ten- 
dency for a simplified color structure in this regime, 
but quantitative estimates remain to be done. Finally, 
we have identified matrix elements closely related to 
two-parton distributions that are suitable for evalua- 
tion in lattice QCD. 

To go from tree level to genuine factorization formu- 
lae, one must be able to sum certain types of collinear 
and soft gluon exchanges into Wilson lines. We argue 
in [12] that for double-scattering processes producing 
color-singlet particles this can be achieved using the 
methods that have been successfully applied to single 
Drell-Yan production [H [I [lOl [11] . This also permits 
the resummation of Sudakov logarithms, with impor- 
tant results sketched in Section [M] For multiple jet 
production the situation is unfortunately more compli- 
cated, as serious obstacles to formulating transverse- 
momentum dependent factorization have been identi- 
fied even for single hard scattering [2^. The produc- 
tion of two electroweak gauge bosons thus emerges as 
a channel where the current perspectives for develop- 
ing the theory look best, and where different aspects 
of multiple interactions can hopefully be explored ex- 
perimentally at LHC. For phenomenological estimates 
of WW production we refer to [64| |65] . 

The cross section for double hard scattering involves 
an integral over the transverse distance y between the 
two partons emerging from each hadron. This includes 
the region of small y, where the two partons can origi- 
nate from the perturbative splitting of a single parton. 
We find that this mechanism, which also affects the 
scale evolution of multiparton distributions, leads to 
serious ultraviolet divergences in the cross section and 
to a double counting problem between single and dou- 
ble scattering. To modify the definition of multiparton 
distributions and possibly the cross section formulae is 
a prerequisite for putting the theory on solid ground. 
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